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Abstract 

Time scales of turbulent strain activity, denoted as the strain persistence times of first and second 
order, are obtained from time-dependent expectation values and correlation functions of lagrangian 
rate-of-strain eigenvalues taken in particularly defined statistical ensembles. Our approach is based 
on a gaussian closure approximation for single-point vorticity statistics, devised to yield a simple 
and computationally efficient method for the evaluation of the strain persistence times. We find 
that softly divergent prefactors correct the usual time-scale estimate of standard turbulence 

phenomenology, a fact that turns out to be of fundamental importance in the emergence of turbulent 
intermittency. The theoretical framework is supported by a straightforward analysis of direct 
numerical simulation data. 

PACS numbers: 47.27.Gs 47.27.eb 
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I. INTRODUCTION 



It is a point of reasonable consensus that further progress in the statistical theory of 
turbulence has been hampered in great part due to the fact that one of its phenomenolog- 
ical pillars - the Kolmogorov-Richardson cascade - is actually a longstanding open issue. 
The usual assumption of vortex stretching as the essential mechanism for the local flow 
of turbulent kinetic energy towards smaller scales has been challenged by the visualization 
of multiscale vortical structures in real and numerical experiments (l|-l5|, and the related 
discovery of geometrical statistics phenomena js, 7|. One may expect that significative ad- 
vances in the derivation of the statistical properties of turbulence will follow from a deeper 
understanding of flow instabilities and their role in the production of coherent structures, 
within more elaborate discussions of the coupled dynamics of vorticity and the rate-of-strain 
tensor. 

A fundamental problem in this context is to determine for how long a given fluid element 
is, in its lagrangian evolution, coherently compressed or stretched by the underlying strain 
field. According to common wisdom js], if s is some measure of the strain strength, such 
a "strain persistence time" can be estimated as T{s) ~ However, this expression for 
T(s) is in fact problematic, since the constancy of sT{s) suggests that weak large scale and 
strong small scale rate-of-strain fluctuations would, respectively, (i) break statistical isotropy 
at small scales and (ii) have no role in the production of coherent structures, as vortex tubes. 

CI 

Both of these implications are at variance with experimental and numerical observations [9|. 

Motivated by the above difficulties, our aim in this work is to show that instead of 
a single time scale T{s), the strain activity is more naturally associated to two distinct 
time scales, which will be denoted as the strain persistence times of first and second order. 
It turns out that these time scales contain divergent prefactors which multiply the usual 
1/s estimate of standard phenomenology, a fact that throws some light on the origin of 
turbulent intermittency and is perhaps related to the well-reported confinement of strong 
vorticity fields in small scale coherent structures. 

This paper is organized as follows. In the next section we address formal definitions of 
the strain persistence times and discuss, by means of a straightforward closure scheme, their 
relation to single-point vorticity statistics. In Sec. Ill, we verify, taking direct numerical 
simulation (DNS) data into account, that our analytical framework, devised to hold in prin- 
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ciple in the small strain domain, incidentally holds for the whole range of strain strengths. 
In Sec. IV, we comment on our findings and point out directions of further research. 

II. STRAIN PERSISTENCE TIMES 

Let Sij = {diVj + djVi)/2 be the (i, j)-component of the lagrangian rate-of-strain tensor. 
Recalling that Sij is traceless due to incompressibility, call the only positive or the only 
negative eigenvalue of Sij{t) by s(t), a piecewise continuous function of time, as indicated in 
Fig. la. 

Independent turbulent flow realizations of s{t) generated, for instance, from some set 
of random initial conditions at t — — oo constitute a large functional space S. Take the 
ensemble C S* of all the profiles s{t) which have s(0) = s for an arbitrarily prescribed 
eigenvalue s. Alternatively, we define the related ensemble A^ of compactly supported func- 
tions s{t) which are identified to s{t) G A^ in the largest neighborhood of t = where s{t) is 
continuous. The functions s{t) vanish out of these neighborhoods. See the sketches in Fig. 
lb. 

In a more formal way, A^ is given as the ensemble of functions s{t) obtained from the 
one-to-one mapping 

As ^ As 

m^sit), (2.1) 

given, for positive t, by 

[ if 3t' e [0,t] I s{t')/s < , 

while for negative t, the time interval [0,t] is replaced, in Eq. (12. 2p . by [t,0]. 

The importance of the ensemble A^ relies on the fact that their elements, i.e., the time- 
dependent strain eigenvalues s{t), have all the same postulated "strain strength" s = s(0). 
Our task, therefore, is to investigate their characteristic time scales and to understand how 
they depend on s. Furthermore, the rationale for the introduction of the additional ensemble 
As, besides A^, is that it leads to expectation values which are less contaminated by errors 
than the ones computed through A,;, while retaining the same (or approximately the same) 
relevant time scales of strain activity. 
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FIG. 1: Figs, (a) and (b) are produced from the direct numerical simulation (DNS) data which 
is addressed in Sec. III. In (a) we show a typical profile of the piecewise continuous function s{t). 
In (b) the compactly supported configurations si{t), S2{t) and S3{t) belong all to the ensemble 
As, with s = —13.6 ±0.1, since si(0) = S2(0) = S3(0) = s. The total time span of the DNS is 
At = 2048 X 10"'^ in arbitrary time units. 

We emphasize that the just given definitions of the ensembles and are strong ideal- 
izations of dubious utility if the interest is to compute statistical averages out of experimental 
or numerical data. The essential difficulty here is that one should work, in principle, with a 



large functional space of turbulent flow realizations. However, by evoking ergodicity 10|, a 



statistically equivalent ensemble A^ can be introduced for practical purposes as it follows: 
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(i) Taking a single three-dimensional turbulent flow realization (e.g., the one obtained in 
a direct numerical simulation) recorded in a time interval of length T, pick up a number 
of (hopefully weakly correlated) lagrangian trajectories; 

(ii) An alternative set of time parametrizations of the lagrangian trajectories is imple- 
mented by conventionally setting time t = at M equally spaced time instants along the 
dynamical evolution; 

(iii) Deflne the ensemble of x M strain eigenvalue proflles s{t) obtained from the 
lagrangian trajectories introduced in (i) and time-parametrized according to (ii); 

(iv) Working with some assigned uncertainty 6s in the deflnition of s, introduce the 
ensemble A^ of functions s{t) derived from the proflles deflned in (iii). 

In the course of analysis, we will be particularly interested in the following time integra- 
tions over flrst and second order expectation values taken in A^, 

his) = f dt{s{t))^^ , (2.3) 

his) = r dt f dt'{5sit)5sit'))j,^ , (2.4) 

where 

5sit) = sit) - (s(t))A, . (2.5) 

Considering, for convenience, positive and negative s as separate cases, strain persistence 
times of flrst and second order, Tiis) and T2is), respectively, can be deflned from (12. 3p 
and (12^ . as 

sTi is) for s < , 



and 



, \sT:^is)Y for s > , 
his) = / (2.7) 



[sT^is)] for s < 

The central idea underlying our discussion is that analytical expressions for T^is) and 
^tis) can be derived from a statistical treatment of vorticity fluctuations, which are gov- 



erned by the lagrangian evolution equation 

duj. 



Ill, 



^ = Sijit)uj + eijkdjfk + vd'^uji , (2.8) 
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where fi is the external force and u is the kinematical viscosity. Let n{t) be the unit 
vector defined along the principal direction associated to the only positive or only negative 
eigenvalue s{t) of the lagrangian rate-of-strain tensor (the twofold orientation ambiguity of 
n{t) is arbitrarily resolved). One gets, after some simple algebra on Eq. fl2.8p . 



dt 



dn 



2s{t)uj + 2ujoj ■ — + eijkUJidjfk + fOJid' 



(2.9) 



where = uj(t) ■ h(t). It is interesting to take the expectation value of Eq. (12. 9p condi- 
tioned to a given time-dependent profile s{t) E Ag. In other words, we change our focus to 
the alternative evolution equation 



dt 



dn 



2s{t){uj )[s] + 2{uu ■ — + eijk{uJidjfk)[s] + v{uJid uji)[s] 



(2.10) 



where ((...))[s] is the self-evident notation for the procedure of conditional averaging. It is 
clear that Eq. (I2.10p is not closed. The second term on its right-hand-side, for instance, is 
related in an intrincate way, through dn/dt, to the velocity gradient tensor and the pressure 
hessian. Neglecting small alignment effects |6i], we take, as a first approximation, that u{t) 
and h{t) are completely uncorrelated in the ensemble of flow realizations which share the 
same arbitrary profile s{t). Resorting furthermore to isotropy, it follows, thus, that 



. 



(2.11) 



The last two terms in Eq. f l2.10p are just the rates of enstrophy injection and dissipation, 
so that their combined contribution is assumed to vanish. Taking Eq. (12. lip into account, 
we are led, after a straightforward integration of Eq. (I2.10p . to 



{uj^{t))is] = (a;^(T))[,]exp 



dt's{t') 



Eq. (I2.12P is now averaged over the configurations s{t) G A 



expressed in path-integral language 



(2.12) 

This can be rethorically 



12j as 



D[sitMm]{^\T))[s^exp 



dt's{t' 



(2.13) 



where D[s(t)]p[s(t)] is the probability measure defined on A^. Considering T to be large 
negative enough (typically much larger, in absolute value, than the large eddy turnover 
time), the squared vorticity expectation value which appears in the integrand of fl2.13p 
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becomes, supposedly, independent on the elements of and on the particular time T. We 
get, therefore, 



(a;^(0))A, =c^ lim /_ D[s(t)]p[s(t)] exp 2/ dt's{t' 



(2.14) 



where = {u'^{T))x^. 

The cumulant expansion method 



13| can, at this point, be straightforwardly applied to 



Eq. fl2.14p to yield, in a gaussian closure approximation (i.e., up to second order in the 
strain fluctuations). 



(? exp 



2/° dt{s{t))x,+2 r dt f dt\8s{t)5s{t'))x^ 

00 J—oo J—oo 



(2.15) 



Noticing that the functions s{t) G are essentially extensions of the corresponding s{t) G A^ 
beyond strain persistence time scales, we can replace, as an approximation, s{t) by s{t) and 
As by As in the right hand side of Eq. fl2.15p . We obtain, thus, from Eqs. (12. 3 p to (12. 7p and 

\ exp l2sTi+(s) + 2[sT2+(s)]2| for s > , 
[ c2 exp {2sTf (s) + 2[sT2"(s)]2| for s < . 
It is convenient to write, up to first order in s, the standard deviation of w(0) in the ensemble 

As as 

a(l + b+s) for s > , 

^ ^ (2.17) 
a(l — 6_s) for s < , 

where a > 0, 6+ and 6_ are arbitrary coefficients and we have used that u and n are 
independent random variables in As, so that (a;(0))A^ = 0. As reported in the next section, 
we have found convincing numerical evidence that there are no higher order corrections to 
Eq. (I2.17p . This is a remarkable fact, still in need of understanding. 

As stated (in rephrased form) in the introductory section, we do not expect vanishing 
strain to have any effect on the statistics of vorticity. Requiring, therefore, that 



lim[sTi±(s)] = lim[sT2=^(s)] = 



(2.18) 



it follows from Eqs. fl^TTB]) and f ETTp that c = a. 

Taking into account, now, the time- reversal symmetry of the fluid dynamic equations in 
the absence of forcing and dissipation terms, a meaningful approximation for the description 
of inertial range processes, we assume that Ti{s) oc Tf (— s) and T2{s) cx T2{—s). We point 



out that this argument is not inconsistent at all with the dissipation anomaly postulated by 



the "zeroth law" of turbulence, that is, the fact that energy dissipation rate per unit volume 



is finite in the inviscid limit z/ — )■ 17(1 . The situation here is analogous to the issue 



on the coexistence of the second law of thermodynamics with microscopic reversibility in the 
statistical mechanics context. Introducing a pair of even functions of s, Fi{s) and -^2(5), and 
proportionality constants g and g', we may write, thus, without loss of further generality, 
that 



T+(s) = -lnFi(s) , Tr(s) = -^lnFi(.) , 

s s 



(2.19) 
(2.20) 



Considering, below, that s > 0, note that 

h{-s) f_^dt{s{t))A_ 



and 



his) 



his) /°^rft(.(t))A. T+(.) 

I' 00 dt dt'{6s{t)6s{t'))^_^ \T,-{-s) 
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T+(.) 



(2.21) 



(2.22) 



his) j'^dtf_^dt'{Ss{t)Ssit'))^^ " 
From Eqs. fl2.16p . f l2.19p and f l2.20p and the even parity of -Fi(s) and F2{s), it turns out that 
Eq. f l2Tfl) holds if and only if 



F^{s) ■ F2{s) = 1 + b+\s\ , 
l^ = l + b \s\ 



(2.23) 
(2.24) 



Substituting the solutions of Eqs. (12:23|) and (^M) for Fi(s) and ^2(5) into Eqs. flCTj) and 
(Km . we find 

1 



THs) 



T+(s) 



9 + 9 



-T,{s) , Tf (s) 



9 

9 + 9 



1 



where 



and 



9 + 9 



Ti(s) = — In 



7^2(5) , T,'{s) 



9' 



9 + 9 



Ms) , 



l + h-\s\ 



T2{s) = —^\n[{l + h+\s\)9{l + h4s\ 



(2.25) 
(2.26) 

(2.27) 
(2.28) 
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It is interesting to remark that once Ti (s) is positive definite, it is necessary to have, 
according to (12.271) . g' > 1 and g'bj^ > 6_ pL8]. This is a simple and well-defined prediction 
from the present formalism. 

In the weak strain regime, one expects that the strain persistence times saturate to the 
large eddy turnover time Tq = ^ L"^ /e, where L and e are, respectively, the typical large 
length scale and the energy dissipation rate parameters of the turbulent fiow. We get, thus, 
from Eqs. (12:251) to ^OTM . 

limT^is) oc g'h+ - 6_ oc Tq , (2.29) 
limT2=^(s) oc gh+ + oc Tq . (2.30) 

Both parameters 6+ and 6_ are, therefore, likely to be proportional to the large eddy turnover 
time Tq. An interesting problem, not touched here, is to find the Reynolds number depen- 
dence, if any, of the dimensionless parameters g, g', b+/To, and b-/To. 



III. ANALYSIS OF DNS DATA 



We have computed statistical averages of fiuid dynamic observables with the help of the 
direct numerical simulation (DNS) database available from the turbulence research group at 
Johns Hopkins University [19!] • An homogeneous and isotropic turbulent fiow with Taylor- 
based Reynolds number R\ ^ 433 is simulated in a periodic cube of linear dimension L = 
271 modeled as a grid of 1024^ lattice points. Viscosity and time step parameters are, 
respectively, u = 1.85 x 10^'* and At = 2 x 10^'^ (the complete simulation record corresponds 
to around one large eddy turnover time. To = 2^° x 10 x At). Further simulation details can 



be found in Refs. 
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2l|. 



Statistical samples were produced in two different ways, according to the particular ex- 
pectation values we were interested to evaluate: 

(i) In order to compute ^^^(s), vorticity vectors and rate-of-strain tensors were defined 
from the velocity gradients taken at grid points 27r(i, 8j, 8A;)/1024 where < i < 1024 and 
^ j,k < 128 are integer numbers, for 10^ frames of equally time-spaced fiow configurations. 
Bins of variable sizes were considered for the sets of positive and negative rate-of-strain 
eigenvalues s. While essentially conventional, our particular bin size choice, 6s = 0.5, 
proved to yield robust results at economical computational costs. 



(ii) For the evaluation of Ii{s) and his) we have considered 10^ lagrangian trajectories, 
each one consisting of 2^° time steps (which are separated in time by 2 x 10^^ arbitrary 
time units). The initial points of the lagrangian trajectories are given by 27r(i, j, k)/10, with 
^ ^ i-i ilk < 10, and have their time evolution determined through a fourth-order Runge- 
Kutta algorithm. The ensemble A^, with uncertainty 5s = 0.1, was generated through the 
procedure previously discussed in Sec. II. 
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FIG. 2: (T^(s) = ^J (a;2(0))^^, i.e., the standard deviation of the projected lagrangian vorticity a;(0), 
is plotted as a function of standardized rate-of-strain eigenvalues s/a±. The standard deviations 
of the positive and negative rate-of-strain eigenvalues s are, respectively, C7+ = 4.63 and (T_ = 6.52. 



In both of the above cases (i) and (ii), the eigenvalues and principal directions of the 
rate-of-strain tensor were computed through an efficient hybrid algorithm which combines 
direct analytical evaluation and the so-called QL algorithm 22 1. 

As it is clear from Fig. 2, the conditional expectation value o"aj(s) is precisely described 
by Eq. f l2.17p . with no additional corrections. The slope parameters 6+ and 6_ associated 
to the right and left branches of ^^^(s), respectively, are 6+ = a^^da^{s) / ds\s>Q = 0.26 and 
b- = a^^dai^{s)/ds\s<o = 0.08. 

In Fig. 3, we show that the expectation value ratios /i(— s)//i(s) and l2{—s) / I-^^s) are 
approximately constant for an extended range of s- values, in agreement with Eqs. f l2.2ip and 
f l2.22p . The measured values of 6+, 6_, g = |/i(— s)//i(s)| and g' = J2(— s)//2(s) are then 
substituted in the logarithmic corrections written down in Eqs. f l2.27p and f l2.28p . which are 
strikingly confirmed from the plots shown in Fig. 4. 
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FIG. 3: The horizontal dotted Unes indicate the mean values —3.84 (= —g) and 15.81 (= g') of 
the ratios Ii{—s)/Ii{s) and l2{ — s)/l2{s), which are, in fact, approximately constant in the range 
3.4 < s < 21.0, the region between vertical dashed lines, where the ensembles are large enough 
for the evaluation of reasonable statistical averages. 



The absolute ratio s)//i(s)| = Ti(— s)/Ti(s) ~ 3.84 is actually expected to be a 
number larger than unity, since, as it is well-known, s is most of the time negative, in other 
words, Ti{—s) > Ti{s). We have in fact verified that the domain of negative s in physical 
space constitutes around 75% of the total fluid volume, which amounts to say that the 
intermediate rate-of-strain eigenvalue is positively skewed {gI. l23|. 
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|s|T2(s) 

FIG. 4: Plots which provide support for the vahdity of Eqs. (j2.6p and (|2.7p with the use of the 
strain persistence times T^{s) and T^{s), as defined from Eqs. (|2.25p - (|2.28p . The dotted hues in 
the plots for have unit slope, while the ones in the plots for /2('S) have slope 2. The upper 

and lower plots in each pair of plots correspond to negative and positive s, respectively. These 
results strongly support the theoretical arguments put forward in Sec. II. The vertical dashed lines 
indicate the range 3.4 < s < 21.0 (see Fig. 3). 

We note, furthermore, that the estimate g' ~ 15.81 is compatible with the relation 
g'b+ > b-, predicted at the end of Sec. II. 
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IV. CONCLUSIONS 



We have devised, within an appropriate statistical framework, a systematic procedure 
for the evaluation of the time scales associated to turbulent lagrangian strain activity. Two 
time scales, the strain persistence times Ti{s) and ^2(3), given in Eqs. f l2.27p and f l2.28p . are 
necessary to accurately reproduce o"a;(s) = y^(w^(0)j^, as determined from DNS data. The 
strain persistence times Ti{s) and T2{s) are introduced as first and second order contributions 
of a gaussian closure approximation. 

While |s|Ti(s) and |s|T2(s) vanish at s = 0, they are slowly divergent at asymptotically 
large |s|, a fact that happens to be a crucial ingredient in the derivation of the linear profiles 
of cr^i^s), clearly indicated in Fig. 2. The divergences of |s|Ti(s) and |s|Ti(s) as |s| — )■ 00 
are actually welcome, since from a phenomenological point of view they imply that the 
coarse grained enstrophy density should grow at smaller scales, where strain fluctuations 
are stronger. Also, it is tempting to conjecture that the behavior of the strain persistence 
times at large s is related to phenomena like the transient formation of unstable vortex 
"pancakes", which tend to roll-up and eventually become vortex tubes, one of the main 
suspected physical mechanisms of turbulent intermittency. 

It is intriguing that the gaussian closure we have implemented here works so well, which 
seems to suggest that the cumulant expansion captures the essential physics of lagrangian 
vorticity fluctuations. It is possible that the second-order truncation in the cumulant ex- 
pansion f l2.15p is actually a fine approximation to the full non-perturbative result, thanks to 
the specific definition of the statistical ensembles A^, which would provide a partition of the 
whole functional space of strain profiles s{t) into subspaces of gaussian stochastic processes. 

We highlight that from a purely theoretical perspective, no considerations have been 
advanced to establish the form of the strain persistence times beyond the first order in s. 
However, we have found that the empirical evaluation of (Ti^{s) does not bring any further 
non-linear corrections into scene, a result which is far from trivial (note that statistical 
isotropy just implies that (cJ^) = 2(sfj), which looks like a necessary but by no means a 
sufficient condition for the specific observed profile of <J^{s)). 

Besides further improvements in the statistical evaluations carried out so far, a natural 
direction of research is to check to what extent the premises and results proposed in this 
work can match the phenomenology implied by promising effective lagrangian simulations 
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of the velocity gradient tensor, as, for instance, the ones derived within the Recent Fluid 
Deformation Closure model 241 125|. 



ACKNOWLEDGMENTS 

The authors acknowledge interesting comments of Arkady Tsinober in the initial stage 
of this work, and important critical remarks by Laurent Chevillard, Bruno Eckhardt and 
Charles Meneveau. This work has been partially supported by CNPq and FAPERJ. 



[1] Z.-S. She, E. Jackson and S.A. Orszag, Nature 344, 226 (1990). 

[2] S. Douady, Y. Couder and M.E. Brachet, Phys. Rev. Lett. 67, 983 (1991). 

[3] D. Bonn, Y. Couder, P.H.J, van Dam and S. Douady, Phys. Rev. E 47, R28 (1993). 

[4] A. Vincent and M. Meneguzzi, J. Fluid Mech. 258, 245 (1994). 

[5] T. Ishihara, T. Gotoh and Y. Kaneda, Ann. Rev. Fluid Mech. 41, 165 (2009). 

[6] Wm.T. Ashurst, A.R. Kerstein and C.H. Gibson, Phys. Fluids 30, 2343 (1987). 

[7] A. Tsinober, An Informal Introduction to Turbulence, Kluwer Academic Publishers (2001). 

[8] H. Tennekes and J.L. Lumley, A First Course in Turbulence, MIT Press (1972). 

[9] To be fair to the usual arguments, the undesirable divergence of T{s) ~ l/l-sj as s — )• has 
not been a source of much concern, since the main focus is, in general, the role of strain in 
inertial range scales. 
[10] B. Galanti and A. Tsinober, Phys. Lett. A 330, 173 (2004). 
[11] S.B. Pope, Turbulent Flows, Cambridge University Press (2000). 

[12] R.P. Feynman and A.R. Hibbs, Quantum Mechanics and Path Integrals, McGraw-Hill Book 
Company (1965). 

[13] R.J. Creswick, H.A. Farach and CP. Poole Jr., Introduction to Renormalization Group Meth- 
ods in Physics, John Wiley & Sons (1992). 

[14] U. Frisch, Turbulence: The Legacy of A.N. Kolmogorov, Cambridge University Press (1995). 

[15] P.A. Davidson, Turbulence, an Introduction for Scientists and Engineers, Oxford University 
Press (2004). 

[16] G. L. Eyink, Nonlinearity 16, 137 (2003). 

[17] Y. Kaneda, T. Ishihara, M. Yokokawa, K. Itakura and A. Uno, Phys. Fluids 15, L21 (2003). 



14 



[18] Straightforward proof of the positive-definiteness of (s): taking a look at Eqs. (j2.3p and 
()2.6p . T^{s) > follows from the fact that {s{t))\^ has always the same sign as the rate-of- 
strain eigenvalue s(0) = s. 
[19] Johns Hopkins University Turbulence Database Cluster, |http://turbulence.pha.jhu.edu/ 
[20] E. Perlman, R. Burns, Y. Li, and C. Meneveau, Supercomputing SC07, ACM, IEEE (2007) 
[21] Y. Li, E. Perlman, M. Wan, Y. Yang, R. Burns, C. Meneveau, R. Burns, S. Chen, A. Szalay 

and G. Eyink, J. Turb. 9, 31 (2008). 
[22] J. Kopp, Int. J. Mod. Phys. C 19, 523 (2008). 

[23] A. Tsinober, E. Kit and T. Dracos, J. Fluid Mech. 242, 169 (1992). 

[24] L. Chevillard and C. Meneveau, Phys. Rev. Lett. 97, 174501 (2006). 

[25] L. Chevillard, C. Meneveau, L. Biferale and F. Toschi, Phys. Fluids 20, 101504 (2008). 



15 



